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Abstract 

We consider possibly observable effects of asymmetric dark matter (ADM) in neutron stars. Since 
dark matter does not self-annihilate in the ADM scenario, dark matter accumulates in neutron 
stars, eventually reaching the Chandrasekhar limit and forming a black hole. We focus on the case 
of scalar ADM, where the constraints from Bose-Einstein condensation and subsequent black hole 
formation are most severe due to the absence of Fermi degeneracy pressure. We also note that 
in some portions of this constrained parameter space, nontrivial effects from Hawking radiation 
can modify our limits. We find that for scalar ADM with mass between 5 MeV and 13 GeV, the 
constraint from nearby neutron stars on the scattering cross section with neutrons ranges from 
cr„ < 10"^^ cm2 to IQ-^^cm^. 
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I. INTRODUCTION 



The characteristics of dark matter (DM) and the nature of its production mechanism 
have so far eluded description. While cosmological observations provide compelling evidence 
for its existence, its mass and the nature of its interactions with the standard model remain 
unknown. One popular hypothesis holds that the DM interacts with standard model particles 
via the weak interaction and is also self-annihilating. In this weakly interacting massive 
particle (WIMP) scenario the correct relic density of DM is a natural consequence of the 
thermal history of the early Universe. 

Alternatively, DM may carry a conserved charge, analogous to baryon number. This 
asymmetric DM (ADM) scenario is motivated by the fact that the DM and baryon densities 
are of the same order of magnitude. The earliest models attempting to relate the DM to the 
baryon asymmetry made use of electroweak sphalerons ^ or out-of-equilibrium decay 
The former often run into tight constraints from LEP measurements. In contrast, the modern 
incarnation of ADM makes use of higher dimension operators to transfer the asymmetry in 
a robust way that is relatively free of electroweak constraints jj]. ADM models prefer to 
have DM mass around a few Ge V (see for example 0-0] ) , which is consistent with hints from 
recent direct detection experiments js], lol. I t is also possible for ADM to have weak scale 
mass , or mass well below a GeV Il2|. 

The DM mass and its scattering cross section with nuclei have been constrained by various 
underground direct detection experiments 13|-[l5|, as well as by particle colliders 16l-ll9|. In 



this article, we study the properties of ADM through its impacts on stellar systems. It has 
long been appreciated that a finite DM-nucleon cross section would result in DM capture 



in stars 



a a ppi 
2JI22 



2Ch22|. In the WIMP scenario, DM annihilation can generate an additional heat 



of white dwarfs 25|, |26| and neutron stars 



source, which may affect stellar formation 23|and evolution 2J], or cause anomalous heating 



27-29 



In the ADM case, DM particles do not 
annihilate and hence provide no additional power for stars. However, since there is no 
annihilation to deplete ADM particles, stars can accumulate far more ADM particles than 
usual WIMPs, which can lead to different effects. For example, it could have implications 
for solar physics 30|-|32| . or change the mass-radius relation of neutron stars 33l-l35|. The 



most extreme possibility is that captured particles can become self-gravitating, forming a 



black hole that will eventually destroy the host stars 36N38 



Recently, constraints on fermionic ADM through the survival of compact stars have been 
discussed in 29, 39|. In a certain class of ADM models, the DM candidate is a boson joi 0]. 
In this paper, we study constraints on this scalar ADM from compact stars. Scalar DM 
particles differ from fermions by spin statistics, which has a significant impact on black 
hole formation conditions. Black hole formation occurs only when the total number of self- 
gravitating DM particles is larger than the Chandrasekhar limit 40|. Fermions obey the 



Pauli exclusion principle, and the Chandrasekhar limit is set by the balance between gravity 
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and the Fermi pressure, while scalar particles have no Fermi pressure to hinder gravity. In 
this case, the lower limit for gravitational collapse is determined by the balance between 
gravity and the pressure induced by the zero point energy, which is much smaller than the 
Pauli pressure experienced by fermions. Therefore, we derive much stronger constraints 
on scalar ADM from compact stars. Since neutron stars have much higher matter density 
and escape velocity than any other stars, we will mainly focus on neutron stars. Note that 



neutron star constraints on scalar DM has been discussed in [36[. In this early work, the 
cross section for the DM capture is set by the geometric cross section, and the DM mass 
range is between 1 GeV and 100 TeV. In this work, we treat the DM-neutron cross section 
as a free parameter and use neutron stars to constrain it. Other important considerations 
that we treat here are the effect of Bose-Einstein condensation, which significantly alters 
the constraint derived on the DM-neutron scattering cross-section, and Hawking radiation, 
which modifies our constraints for high mass DM. We also explore a wider DM mass range 
and use the most recent neutron star data. 

This paper is organized as follows. In Sec. II, we discuss the Chandrasekhar limit for 
fermions and bosons. In Sec. Ill, we discuss ADM capture in neutron stars. In Sec. IV, we 
discuss thermalization, condensation and black hole formation of captured scalar ADM in 
neutron stars. In Sec. V, we discuss Hawking radiation and destruction of the host neutron 
star. In Sec. VI, we discuss observational constraints. We present our conclusions in Sec. 
VII. 



II. CHANDRASEKHAR LIMIT 



First, we review the derivation of the Chandrasekhar limit for a system of fermions. 
Suppose there are N fermions of mass m distributed in a sphere with radius R, so that the 
number density of fermions is approximately N/R^. Because of the Pauli exclusion principle, 
the average distance between two fermions is ~ R/N^^^. The uncertainty principle requires 
that each fermion have Fermi momentum p ~ N^^^/R. If the total number is small and 
m > p N'^^'^ / R, the system is in the nonrelativistic limit. The average energy per fermion 
is 

GNm^ 1 /A^V3~ 
E ^ + - , 1 



R m\ R 

where G is Newton's constant. Once the gravitational and Fermi pressures reach equilibrium, 
the system can have a stable spherical configuration with radius 

R (2) 
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As increases, the radius shrinks and the Fermi momentum increases; eventually fermions 
become relativistic with total energy 

GNm'^ 

E + . 3 

R R ^ ' 

If the total number of the fermions increases beyond the limit 

where Mpi = 1.2211 x 10^^ GeV is the Planck scale, the gravitational energy will dominate the 
total par ticle energy and gravitational collapse will occur. This is the famous Chandrasekhar 



limit 40 



Now we discuss bosons. Similar to the fermion case, the gravitational collapse occurs when 
particles are relativistic. But the bosonic system is significantly different from the fermionic 
system because it has no Fermi pressure to hinder gravity. Since the bosons are confined 
inside a sphere with radius R, they have zero point energy 1/R due to the uncertainty 
principle in the relativistic limit. Therefore, the typical energy for a boson in a sphere of 
radius R is 

GNm'^ 1 

Again, the radius cancels in the critical limit. In this case, the Chandrasekhar limit is 

\ m J \ m J 

Comparing Eq. (jl]) and Eq. ([H]), we can see that for a given particle mass, a particle that 
obeys Bose-Einstein statistics will experience gravitational collapse much more readily than 
a particle that obeys Fermi-Dirac statistics. 

When the total number of DM particles accumulated in a neutron star surpasses the 
Chandrasekhar limit, the captured DM particles collapse to a black hole and destroy the 
host neutron star. Therefore, observations of old neutron stars can be used to constrain 
the DM-neutron scattering cross section. Since bosons have much smaller Chandrasekhar 
limit than fermions, we can obtain stronger limits on bosonic DM. In this work, we take 
typical neutron star parameters M„ = 1.44 Mq, i?„ = 10.6 km and the central density 



Pb = 1.4 X 10^^ g/cm^ [29 



41 



III. CAPTURE OF ASYMMETRIC DARK MATTER IN NEUTRON STARS 



The accretion of DM onto stars has been studied in 20-22 . In this section, we review 



the basic formulas for the capture of asymmetric DM in neutron stars. 
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In the absence of annihilation, the number of DM particles in a star is determined by the 
differential equation 

Cb, (7) 



dNx 



dt 

where Nx is the total number of DM particles in the star and is the DM capture rate 
through scattering with baryons. In this study, we assume there is no symmetric component 
and neglect DM annihilation completely. In some ADM models, a portion of the symmetric 
component can be regenerated by oscillation effects For those models, bounds 

derived from the survival of the neutron stars can be weaker. 

Additionally, we ignore the self-capture effect. Dark matter may have a sizable self- 



interaction that leads to self-capture [42] . However, the self-capture saturates when the sum 
of the individual self-scattering cross sections becomes larger than the geometrical area over 
which the DM particles thermally distribute. As we will show in the next section, due to 
the large baryon density the captured DM particles are thermally distributed in a very small 
region of radius ~ 1 m in the core of the neutron star. We have checked that the baryonic 
capture always dominates the DM accretion process in neutron stars for the parameter space 
of interest, even if we take the upper limit of the self-scattering cross section allowed by the 



eUiptical halo shape bound [43 . 



A. Capture rate 



The number of DM particles in the neutron star is set by Cb- Since neutrons are 
degenerate in the neutron stars, capture can occur only when the momentum transfer is 
larger than the difference between the Fermi momentum and the neutron momentum. As 
we will show, this will affect the capture efficiency significantly for DM with mass less than 
1 GeV. For larger DM mass the effect is negligible, because the momentum transfer is always 
sizable. 



The accretion rate is given by 21 

Cb = 47r 



^" ,dCB{r 



dV 



where i?„ is the radius of the neutron star and the capture rate per unit volume for an 
observer at rest with respect to the DM distribution is given by 



dCBir) 
dV 



6 f(r)^ 
nx{r)nB{r)^——{v(TxB) 



1 - 



1 - exp (-5^ 
52 



(9) 



Here nx{r) is the ambient DM number density; nsir) is the number density of the stellar 
baryons; v is the DM velocity dispersion around the neutron star; v{r) is the escape velocity 
of the neutron star at the given radius r; axB is the effective scattering cross section 
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between DM particles and nucleons in the neutron star; and ^ takes into account the neutron 
degeneracy effect on the capture. The factor is given by 



where = mx/rriB and fi^ = {fi — l)/2. 



3 v{rY fi 
2 



2 ' 



(10) 



Now we specify the factor C,- AH energy levels below the Fermi momentum pp have been 
occupied. During the scattering process, if the momentum transfer to the neutron is larger 
than pf, the scattered neutron can be excited above the Fermi surface. In this case, all 
neutrons can participate the capture process, and the capture efficiency is ^ = 1. On the 
other hand, if the momentum transfer 6p is less than pp, only those neutrons with momentum 
larger than pp — 6p can participate in the capture process. The fraction of these neutrons 
is ^ 6p/pF, so we can approximate a.s — Sp/pp. Depending upon the momentum transfer 
6p, we can parameterize ^ as 



^ = Min 



5p_ ■ 

Pf 



:iv. 



When the DM particle approaches the surface of the neutron star, its velocity is close to 
the escape velocity. Hence, the typical momentum transfer is Sp ~ ^J^rrirVesc-, where = 
'^x'^slij^x + ^^b) is the reduced mass, and typically v^sc — 1-8 x 10^ km/s. The Fermi 
momentum \spp ^ (Svr^pB/ms)^'''^ — 0.575 GeV for ps = 1.4 x 10^^ g/cm'^. Therefore, ^ — 1 
for all mx ^ 1 GeV. In contrast, the capture rate is suppressed by a factor ~ mxVesc/PF if 
the DM mass smaller than the neutron mass. 

To estimate the capture rate, we take the conservative limit that v{r) = v{Rn) = Vesci 
and we assume that nx{r) and nB{r) are independent of the radius; thus, the total capture 
rate can be simplified to 



Cb ^ \f^^^{vaxB)^NB 

V TTTUx f 



exp 



-B' 



52 



(12) 



where A^^ is the total number of neutrons in the host star. When B^ ^ 1, the term in the 
square bracket is close to 1; for typical values Vesc — 1-8 x 10^ km/s and v ~ 220 km/s, and 
a DM mass smaller than 9.4 x 10^ GeV, this condition is obtained. If DM has mass larger 
than ~ 9.4 x 10^ GeV, a lower probability to lose enough kinetic energy to be captured after 
single scatter results 29|] . In our numerical work, we use the full expression of Eq. f[T^ . 



B. Total Number of ADM in Neutron Stars 



Here we only consider DM particles captured by the neutron star itself and neglect those 
the neutron star can inherit from its progenitor. Compared to the neutron star phase, the 
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progenitor usually has much lower density and shorter lifetime which results in lower capture 
efficiency. The total number of DM particles captured by the neutron star is given by the 
solution of Eq. ([7]) 

Nx = CbI. (13) 

To evaluate Cb,^g note that if the sum of individual nucleon-DM scattering cross sections is 
larger than the geometric surface area of the star, the capture rate will saturate. Therefore, 
the capture rate increases with the cumulative nucleon-DM scattering cross section atot = 
NsO'n, where cr^ is the DM-neutron elastic scattering cross section, as long as atot is smaller 
than a geom = ^ri?^; that is, we can constrain the individual scattering cross section cr^ as 
long as (T„ is less than or equal to <Jmax = t^R^INb- Taking typical neutron star parameters 



a_ = 2.1 X 10-^^ cmM — ^ ) ( ^^^^ ) , (14) 



M„ = 1.44 Mq and i?„ = 10.6 km, we estimate the maximum cross section as [29 

2 

10.6 kmy V 

and the effective cross section is given by 

axB = Min [cr„, amax] ■ (15) 

Note that since we consider scattering off only one nucleon, this scattering can be regarded 
as either spin-dependent or spin-independent. 

Now we can estimate the total number of ADM in the neutron star at a given time, using 
generic parameters Vesc = 1-8 x 10^ km/s, v = 220 km/s, and Nb — 1.7 x 10^^. In the regime 
mx ^ 1 GeV, we have — 1, which gives 

/V ~ 2 3 10^4 /100GeV\ ( PX \ ( crxB \ ( t \ 
^ ~ V / V 103 GeV/cm3 ) U.l x lO'^s cm^ ) \ years / ' ^ ' 



When the DM mass is less than ~ 1 GeV, the degeneracy effect on the capture process is 
important so that ^ ~ \/2mxVesclVF-, and we have 

-'^1x10^^/^ \{ \{ ^ \ 

^ ' V GeV/cm3 y V 2- 1x10-45 cm2y' V 10^° years; ^ ' 

It is interesting to note that the DM number does not depend on the DM mass in the second 
case. 

In the above derivation of Nx , we have assumed that the evaporation effect is negligible for 
the DM. Now we estimate the DM mass scale below which the evaporation is relevant. Since 
energy states below the Fermi surface are occupied, only those neutrons with momentum 
above pp can transfer kinetic energy to the DM. Since T <^ pp for the neutron star, the 
number of these free neutrons is order ~ 10~* smaller than that of the neutrons in the Fermi 
sea. So the scattering probability for the DM evaporation is highly suppressed. Furthermore, 
compared to the sun, neutron stars have much higher density and deeper gravitational wells. 
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so it is much more difficult to accelerate trapped DM above the escape velocity through 
interactions with neutrons. To evaporate from the neutron star, the DM has to gain enough 
energy such that its velocity is larger than the escape velocity of the neutron star. Because 
of the degeneracy effect, the typical energy transfer from the free neutron is ~ T; so the 
evaporation effect is relevant only when the DM mass is less than ~ 2T /vl^^ ~ 48 eV for 
T = 10^ K, which is much below the lower mass limit of our constraints ~ 2 keV in the most 
optimistic case. Therefore, we can safely ignore the evaporation process. 



IV. ASYMMETRIC SCALAR DARK MATTER IN NEUTRON STARS 



A. Thermalization 



When DM particles are captured by the neutron star they lose energy via scattering with 
neutrons, and soon attain thermal equilibrium with the star. To estimate the thermalization 
time scale we calculate the DM energy loss rate: 

dE 

= -inBCTnV 5E, (18) 

where ub is the neutron number density in the center of the neutron star, 6E is the energy loss 
of the DM particle during each scattering event, and we use ^ defined as in Eq. f lTTj) to param- 
eterize the neutron degeneracy effect on the DM thermalization process. The typical velocity 
and the momentum transfer 5p = \/2mrV fully determine ^. However, unlike the capture case, 
where the velocity is set by the escape velocity Vesc, the thermal equilibrium of DM particles 
and neutrons now sets v ~ ^J^Ethjmx-, where Et^ ~ 3T/2 is the energy after thermalization. 
In the case of > 1 GeV, 5p ~ ^iubv ~ 2.1 x 10"^ GeV(T/105 K)^/2(100 GeV/mx)^/^ 
which is much smaller than "pp ^ 0.575 GeV. For mx ^ 1 GeV, the momentum transfer is 
given by bp ~ \f2mxv ^ 6.8 x 10^^(T/10^ K)-'^/^(mx/0.1 GeV)^/^, which is again less than 
pp. Therefore, the neutron degeneracy effect reduces the DM thermalization efficiency over 
the entire DM mass range, and ^ is everywhere given by ^ ~ bpjpp. 

To estimate the thermalization time scale, we solve Eq. f lTSj) and get 

. . (19) 

In the limit of vnx ^ 1 GeV, the thermalization time scale can be further simplified to 
tth — "\/2^xPF/(12m^r;,BO"„T). Taking typical values, we see 



/ mx \2 /2.1 X 10-^5 cm2\ /lO^ K\ , , 
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If DM mass is less than 1 GeV, the thermahzation time scale is given by 

- X 10- years (— — j ( j (— j . (21) 

To derive constraints on scalar ADM from black hole formation, we will assume the captured 
scalar ADM follows the thermal distribution in the neutron star. This is only true when tth 
is less than the neutron star age ~ 10^° years. As we can see from Eqs. (120!) and (12T!) . light 
DM easily satisfies this condition. For heavy DM, t^h is not always less than the neutron 
star age. In the following discussion we first assume the DM reaches thermal equilibrium 
with neutrons, and then we check the consistency of this assumption. 

After attaining thermal equilibrium, captured DM particles drift to the center of the star 
and form an isothermal distribution with the typical radius 

T 100 GeV y/^ 
10^ K mx J ' ^ ^ 

We can see that the captured DM particles very quickly occupy a very small region near the 
neutron star core. 



B. Self-gravitation and Black Hole Formation 



If the DM density is larger than the baryon density within the thermal radius vth, the DM 
particles can become self-gravitating. For a total DM mass Mx = Nxttlx within a thermal 
radius rth, this condition is 

3Mx 



PB- (23) 



Therefore, the DM becomes self-gravitating once the total number of DM particles is larger 
than a critical number 

100 GeV\ T ^ 



A,„,,.4.8xl0"^^^j . (24) 

Recall the upper limit for the bosonic system given in Eq. (E]) above which the zero point 
energy cannot prevent gravitational collapse 

iV^r" - 1-5 X 10=^^ ( Mi!£X^ ' . (25) 
\ mx J 



Thus, if the scalar ADM thermalizes and the mass satisfies rrix ^ 10^'^ GeV(T/10^ K)^, 
we always have Ngeif ^ ^chT'- ■'■^ ^^i^ case, gravitational collapse occurs as soon as DM 
particles become self-gravitating in neutron stars. 



9 



C. Bose-Einstein Condensation 



In the above discussion, we imphcitly assumed that all captured scalar ADM particles 
followed a Maxwellian velocity distribution. At the extreme densities we are considering 
here, however, this minimal assumption is not necessarily satisfied. In particular, ensembles 
of bosonic particles at high densities exhibit novel statistical properties. If the central 
temperature of the neutron star falls below the critical temperature to form a Bose-Einstein 
condensate (BEG), the particles in the ground state condense and no longer follow the 
thermal distribution. We will now show that for light ADM this condensation increases the 
density and reduces the restriction on self-gravitation to such an extent that the number of 
ADM particles necessary for self-gravitation is less than the bosonic Chandrasekhar limit. 
Thus, gravitational collapse is set by N^^°^. 

To check this sequence of events, we begin by noting that for a given bosonic DM number 
density nx-, the critical temperature to form a BEC is given by 



= 



mx 



nx 



C(3/2), 



2/3 

(26) 



where C, is the Riemann-Zeta function, C(3/2) ~ 2.612, and nx = 3A^x/(47rr^/j). To see how 
likely it is that the captured ADM will form a BEC in the neutron star, we can estimate the 
critical ADM number as 

where we have used Eq. (1221) . Therefore, if the total number of captured ADM in the 
neutron star is larger than 1.0 x 10^^ (T/10^ K) , some of captured ADM particles will go 
to the ground state and form a BEC. This condition can be satisfied for a neutron star with 
relatively low central temperature as indicated by Eqs. ( |T6l) and ( |T71) . 

For T < Tc, the BEC forms and the number of particles in the condensed ground state is 

\ 3/2" 



N^x = Nx 



1 



-Nx- 1.0 X 10^^ ( ) . (28) 



Since these ground-state particles effectively have zero temperature, they sink deep into the 
core of the neutron star. We can estimate the radius of distribution of the ground state by 
requiring the zero point energy equal the gravitational energy 

f ^ V'^ , /100GeV\'/' , , 

^BEC = - 1-5 X 10-^ cm . 29 

\8nGmi^PB J \ nix J 

This is much smaller than r^^, which indicates a much higher DM density. Thus, the ground 
state itself may become self-gravitating. The critical number for the self-gravity of the DM 
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particles in the condensed state is 



^0 _ 47r psrl^c , , ^ ,23 / 100GeV\ 

\ mx J 



= - 1.0 X lO^M ^^^^^ . (30) 



Once the number of DM particles in the ground state is larger than N^eij, these ground-state 
particles become self-gravitating. Since N^eif less than Nseif, the onset of self-gravity is 
marked by N^eif instead of A^sei/ in conditions where a BEG forms. As indicated above, this 
leads to qualitatively different behavior as compared to the case when a BEG does not form, 
since now N^eif can be less than N^^°^. If this is the case, as soon as a condensed ADM 
system reaches the Ghandrasekhar limit it will undergo gravitational collapse. 

For this effect to be important, A^^ has to grow larger than the Ghandrasekhar limit 
for a bosonic system, so that the condition for black hole formation of the BEG becomes 
^ ^cha^- -^y using Eq. fl28|) . we get a lower limit on the total DM number Nx, 



Nbec = N'ctr + 1-0 X 10=^^ 



T 



105 K 

2 



We can see that the value of the right-hand side of Eq. fl3T]) is less than Ng^if if the DM mass 
mx < 1.9 X 10^ GeV (4.7 x 10^ GeV) for a central temperature T = 10^ K (10^ K). In the 
situation where this condensation can occur, the BEG shortens the time scale for the black 
hole formation for the scalar ADM in the low-mass range. 



V. HAWKING RADIATION AND DESTRUCTION OF THE HOST STAR 



During the collapse process, the gravitational contraction releases energy which can be 
absorbed by neutrons through DM-neutron scattering. This cooling mechanism is so efficient 
that eventually the DM sphere collapses to a black hole 36|]. Once a black hole is formed 
at the center of the neutron star, it will rapidly capture the baryonic matter of the neutron 
star. Hawking radiation will also be active, reducing the mass of the black hole and possibly 
heating the remaining DM. Finally, the black hole may also consume the ambient DM 
particles, which can be crucial for the stability of the black hole. Here, we analyze the 
relative contributions of these effects, and we see that in the majority of our parameter 
space the black hole will grow and eventually consume the host neutron star. Once the 
physics of the accretion is properly considered, we find that Hawking radiation could be 
important for high- and intermediate-mass ADM. 

For a black hole with mass Mbh , the differential equation that governs the rate of change 
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of mass is 



dt - ^"^^ J ^^"^ - 15360vrG2Ml^ + J " ^^^^ 

The second term of Eq. (1321) represents the Hawking radiation rate, while the third term is 
the accretion rate of ambient DM particles. The first term of right-hand side of Eq. fl32|) 
is the Bondi-Hoyle accretion rate, in which Vs = \fdPJdp is the sound speed and is the 
accretion eigenvalue for the transonic solution. To determine Vs and A^, we characterize 
the equation of state of the neutrons by P = Kp^ ^ where K and 7 are constant. For a 
nonrelativistic degenerate neutron gas, which is a good approximation for neutrons in the 
neutron star, we have 7 = 5/3 and K = S'^^^n^^^ / {5m^^^) 40|. We estimate the sound speed 
as Vs = \/ K'yp'^^^ ~ 10^ km/s, where we take p ~ 1.4 x 10"^^ g/cm^. The accretion constant 
is given by A, = (l/2)(^+i)/(27-2) [(5 _ 3^)/4]-(5-37)/(27-2) _ o.25 



dMBH , . (GMbhY 1 ( dMBH 



40 



A. Black Hole Mass without Bose-Einstein Condensation 



For large ADM mass, we found above that Ngeif < Nbec and a black hole forms without 
the assistance of a BEC. After formation of the black hole, the neutron star continues 
capturing DM particles. These newly captured DM particles eventually sink to the center of 
the neutron star and distribute themselves within r^/j. In principle, the black hole can increase 
its mass by capturing these additional DM particles. However, we find that this capture 
rate is very small and {dMBH I dt) dm is negligible. This is because, for a nonrelativistic 
particle moving towards the black hole, its impact parameter must be less than hmax = 



AGMbh /voo [4Q|] to penetrate the angular momentum barrier and fall into the black hole. 
Here, foo is the particle's velocity when it is far away from the black hole. 

Taking Mbh ~ ^xNsdf and Voo ~ A/3T/mx, we can estimate hmax as 

, /lO TeVA f T\ 
bmax ~ 1.6 X 10-^ cm h (33) 



mx J \10^ J 

which is much smaller than the thermal radius rth — 2.4 cm(T/10^ K)(10 TeV/m^)^^^. 
Therefore the majority of DM particles captured after formation of the black hole do not fall 
directly into the black hole. The remaining DM particles orbit the black hole at a distance 
of order rth] the black hole gains mass from these particles at a rate set by the coUisionless 



spherical accretion approximation |40|. We find that the DM accretion rate is much less 
than the baryon accretion rate, so we can safely ignore the [dMBu/ dt)£,M term in this case, 
and we obtain a critical initial black hole mass 

M^^l ~ 1.2 X 10^^ GeV. (34) 

Without a BEC, the initial black hole mass is Mbh ~ Nseiffrix- If we demand Nseif^ix ^ 
M'g^, we find that for mx < 2.6 x 10^ GeV (T/10^ K) the Hawking radiation has a 
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longer time scale than the accretion process. Hence in this mass range the black hole will 
continuously accrete baryonic matter until the neutron star is consumed entirely. 



B. Black Hole Mass with Bose-Einstein Condensation 

For low-mass ADM, particles in the BEC ground state form a black hole. We must check 
the mass above which the black hole evaporates. If we naively ignore the term {(IMbh / dt) dm 
and demand Mbh ~ ^xN^^°^ > ^bhi ^"^^ i\iaX the black hole mass increases only for 
the DM mass less than ~ 13 GeV. If so, the constraint is valid for mx ^ 13 GeV. But 
in contrast to the non-BEC case, {(IMbh / dt) dm inay have an important effect on the black 
hole mass evolution, and we find that the bound can be sensitive to masses higher than 
~ 13 GeV. We detail our reasoning below. 

Since the black hole forms only from ADM particles in the ground state, the remaining 
ADM particles follow an isothermal distribution with a radius rth. As discussed above, the 
thermally distributed DM particles do not fall into the black hole, and so the phase space 
of the non-BEC state is still completely occupied. Hence, if any more ADM particles are 
introduced to the thermal region, a new BEC ground state must form in the center of the 
star. In this way, the introduction of more DM particles into the thermal radius essentially 
forces the formation of a BEC ground state. Before and after the mini black hole forms, the 
neutron star continuously captures ADM particles. All of the captured ADM particles will 
eventually thermalize, sink to the center of the neutron star, and prompt the formation of a 
new BEC state. 

If the thermalization time scale is shorter than the evaporation time scale of the mini 
black hole, the black hole can always efficiently accrete ADM particles in the new BEC 
state. Taking the initial black hole mass as Mbh ~ ""^xA^cvia" ~ ^/{Gnix), we can 
conservatively estimate the black hole evaporation time scale in the absence of particle 
accretion as t^aw — l53607rG'^M^jj/3 ~ 5 x 10^ years (100 GeV/mx)^. This is much longer 
than the thermalization time scale given in Eq. fl20l) . Since all newly captured ADM particles 
eventually go to the ground state after the amount of time it takes them to thermalize, the 
rate at which ADM particles fall into the BEC ground state per unit time is given by the 
capture rate Cb- To check whether or not the ADM particles in the BEC state feed the 
black hole efficiently we calculate the maximal impact parameter bmax and compare it with 
the distribution radius tbec- With Mbh ~ ^/{Gmx) and ~ ^/{mxrBEc), we have 

brnax ~ 'irBEC, (35) 

where tbec given by Eq. fl29l) . Since bmax ^ rBEC, we see that the black hole can efficiently 
consume the BEC. This occurs at a rate given by {dMnn / dt)DM ~ rnxCB- 
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This new source of accretion overwhelms the Hawking radiation, which is emitted at the rate 
{(IMbh I dt)j^^ ~ 10^'^ GeV/year (mx/100 GeV)^. Thus, in the BEG case our constraints can 
be sensitive to scalar ADM with mass much higher than the naive estimate, mx ~ 13 GeV. 

In this discussion, we have assumed that the black hole formation and Hawking radiation 
do not destroy the BEG state. However, one may note that the black hole radiates as a 
blackbody at a very high temperature. Thaw = {87^GMbh)~^ ~ mx/Svr. If this Hawking 
radiation is emitted solely as relativistic particles of the standard model such as photons 
and neutrinos, these particles will heat the neutrons which in turn heat the ADM. We may 
estimate the change in neutron temperature by assuming that all of the initial rest mass 
energy of the black hole goes into heat. We see that in this case equipartition of energy 
requires that AT ~ mxN^^^ /Nb — 1-0 x 10^^ K (10 GeV/mx), and there is no change in 
the thermal ADM distribution. 

On the other hand, if the energy produced by the Hawking radiation is efficiently trans- 
ferred to DM particles, the situation may change. For example, the DM may couple to 
some light mediator particles, which can be produced by the Hawking radiation. These light 
mediators may heat the ADM directly and transfer the black hole energy to the thermal 
energy of ADM particles. In this case, the calculation is rather complicated and a detailed 
analysis is beyond the scope of the current paper. Here we give a conservative estimate by as- 
suming that all of the initial black hole mass goes to the thermal energy of DM particles. The 
equipartition of energy gives AT ~ mxN^chTl^th ^ 1.7x 10^^ K (10 GeV/mx) (T/10^ K)"^ 
where the number of remaining non- BEG particles is Nth = 1-0 x lO^^ (T/IO^ K) , from Eq. 
(!3T|) . This heating will greatly affect the thermal distribution of the scalar ADM particles. 
Since rth oc a/T, the phase space will expand greatly, and the newly captured ADM particles 
do not form the BEG state. Thus, the black hole constraints are lifted in this case. 

In practice, for heating processes to occur, ADM particles must couple to light mediator 
particles, since the black hole can only produce particles with mass much less than Thaw ~ 
mx/SiT. On the other hand, the presence of the light mediator may also give rise to DM self- 
interactions, and the observed ellipticity of DM halos places a lower bound on the mediator 



mass [43|. As an example, if mx ~ 13 GeV, the halo shape constraint requires the mediator 
mass be larger than ~ 40 MeV 4J|. Thus, for scalar ADM with mass mx ~ 13 GeV the 
existence of a mediator with mass in the range ~ 40 MeV — 0.5 GeV may help in evading the 
black hole formation constraints, while remaining consistent with the halo shape bound. We 
can see that the constraints for mx ^13 GeV in the BEG case are rather model-dependent, 
and we specify these regions in our plots. 

Now we estimate the destruction time scale. If the initial black hole mass exceeds the 
critical value M^}^, the time scale to destroy the neutron star is set by the Bondi-Hoyle rate 

i - r./\. ^ 2.3 X 10-^ s ( ^\ , (37) 
txG^PbMbh \MhH) 

where we take Vs ~ 10^ km/s and Mq is the mass of sun. Therefore, the characteristic time 
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scale is t ~ 17 years (mx/100 GeVr^ (10^ K/Tf/^ and 5.4 x 10^ years (mx/GeV) for the 
non-BEC case and BEG case, respectively, which we note is much shorter than the typical 
old neutron star age ~ 10^° years. 



VI. OBSERVATIONAL CONSTRAINTS 



So far, we have gone through the conditions for the captured scalar ADM particles to form 
a black hole. Such a black hole can destroy the host neutron star. We can see that these 
conditions are easily satisfied if the ADM particles have a sizable scattering cross section 
with neutrons. However, we observe many old neutron stars near the solar system and in 
globular clusters. Therefore, these observed old neutrons stars constrain the ADM-neutron 
scattering cross section. In this section, we proceed to derive upper bounds on the ADM- 
neutron scattering cross section for a given ADM mass. In Sees. IVI Al (IVI B|) we first derive 
constraints on (T„ by assuming fiducial neutron star parameters for the case without (with) 
Bose-Einstein condensation. In Sec. IVI Ct we apply these bounds to observed neutron stars 
in the solar neighborhood and in the globular cluster M4. 



A. Constraints on DM-neutron cross section without Bose-Einstein condensation 



We first discuss constraints in the absence of BEG formation; this gives rise to conservative 
constraints in the low-mass range but is more stringent for high masses, where BEGs do 
not readily form. Since in this case Nseif > N^°hr ^^i mx < 10^^ GeV(T/10^ K)^, 
gravitational collapse occurs as soon as DM particles start self-gravitating. In order to 
avoid the destruction of neutron stars, we demand Nx < Ngeif and get an upper bound on 
DM-neutron scattering cross section cr„. For DM with mass mx ^ 1 GeV we derive the 
following bound on cr„ using Eqs. f|T6|) and (1241) and requiring that the host neutron star is 
not destroyed over the course of its lifetime t 

.^ 48 2 /^lO^ GeV/cm3 10^° years \ / 100 GeV T V'^ 

an < 4.4 X 10"^^ cm^ '- — . 38 

\ Px t \ mx IQ^'K ^ ' 



For mx < 1 GeV, the bound on cr„ is 



.„<,,xlO-«c,„M^^^^.^^j^^j . (39) 



5/2 / rp \ 3/2 



Since the capture rate saturates when the DM-neutron scattering cross section is larger than 
amax — 2.1 X 10~^^ cm^, the upper bound on cr^ is only valid when the value of the right-hand 
side of Eq. (l38l) and Eq. (|39l) is smaller than amax- Because the capture and thermalization 
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processes do not distinguish between spin-dependent and spin-independent cross sections, 
the bound on (T„ apphes to both cases. 

We depict the constraints from requiring neutron star survival in the left panel of Fig. ([1]) 
with various values of the DM density. We take the central temperature as 10^ K and neutron 
star age to be 10^° years. The most prominent qualitative feature, the sharp vertical cutoff, 
corresponds to cross sections cr„ ~ 2.1 x 10~^^ cm^. Here, the geometric cross section limits 
the capture of DM particles, so we cannot constrain the interaction cross section for the mass 
below the cutoff. Furthermore, without a BEC, we can constrain scalar ADM with mass 
mx ^ 100 MeV only if the DM density px ^ 10® GeV/cm^. One may find regions with such 
high DM density near the galactic center. As we discussed before, if the DM mass is less than 
~ 1 GeV, the capture rate is reduced due to the neutron degeneracy effect, which is indicated 
by the change in slope of the curve with px = 10® GeV/cm^ when mx ^ 1 GeV. For DM 
mass mx ^10® GeV, the factor in Eq. (fT2!) falls below 1 and the capture probability 
through a single scatter becomes lower, as shown by the bump in the left panel of the figure. 
In the hatched region, DM particles cannot thermalize with the surrounding neutrons within 
the age of neutron star as discussed in Sec. IV A. In this region, the captured DM does not 
necessarily distribute within the small thermal radius of the neutron star core, and the bound 
does not apply. In the cross-hatched region, the initial black hole mass is so small that it can 
evaporate due to the Hawking radiation as shown in Sec. V A. In the square-hatched region, 
the Hawking radiation may interfere with DM accretion, and the black hole will evaporate 
in some case as discussed in Sec. V B. 

We can also estimate the analogous bound from white dwarf stars. Typical white dwarf 



parameters are [45|: mass ~ 0.7 Mq, radius ~ 6.3 x 10'^ km, density ~ 10® g/cm^, central 
temperature ~ 10^ K and escape velocity ~ 6 x 10^ km/s. White dwarfs are composed of 
carbon and oxygen; we make the conservative assumption that the white dwarf is entirely 
composed of carbon. We find that the white dwarf bound on (T„ is about 9 orders of 
magnitude weaker than the limit derived for typical neutron stars. 



B. Constraints on DM-neutron cross section with Bose-Einstein condensation 



If the captured scalar ADM does form a BEC, the constraints on the DM-neutron cross 
section become stronger. This is because Nbec ^ A^se«/ for the DM mass less than a few 
TeV, dependent on the central temperature T. In this case, self-gravitation and gravitational 
collapse will occur more quickly due to the heightened density of the ADM. 
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(GeV) (GeV) 

Figure 1. Regions (gray) of DM- neutron scattering cross section in which accumulated scalar ADM 
forms a black hole. We have px = 0.1 GeV/cm^ along the dotted red lines, which is approximately 
the DM density in the solar neighborhood. In the left panel, we assume a BEC does not form and 
all captured DM particles become self-gravitating and collapse. In the right panel, we assume a 
BEC does form and DM particles in the BEC ground-state collapse and form a black hole. We 
take the neutron star age to be 10^'' years with a central temperature 10^ K. In the diagonally 
shaded regions, DM particles cannot thermalize with neutrons within the age of neutron star. In 
the diagonally cross-hatched region, the black hole can evaporate due to the Hawking radiation. 
In the square-hatched regions. Hawking radiation may interfere with DM accretion by heating the 
thermalized DM, and in some cases the black hole can evaporate (see discussion in Sec. V). In these 
hatched regions, the bounds are lifted. The black regions are excluded by recent CDMS results 
(spin-independent) [3]. 



For mx ^ 1 GeV, the requirement Nx ^ Nbec gives 



(T„ < 9.1 X 10-^^ cm^ 
1.5 X IQ-^ 



10^ GeV/cm^ 10^° years ^ f nix 

Px t 
100GeV\^ / T 



100 GeV 



mx 



105 K 



(40) 
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In the case of mx ^ 1 GeV, the upper hmit on (T„ is given by 

lO^GeV/cm^ 10^° years 



(T„ < 6.2 X 10"^^ cm^ 



Px t 



X 



1.5 X lOM ^i^V + ^ ^ ^ 



mx 



] / ^ Y 



(41) 



Again, we derive bounds by requiring that the host star is not destroyed by the gravitational 
collapse of the ADM. 

In the right panel of Fig. ([1]), we display the DM-neutron scattering cross section with 
various values of the DM density that satisfy Nx ^ Nbec- For mx < 10 GeV, the BEG 
forms before the Ghandrasekhar mass is reached, while for mx ^ 10 GeV, the Ghandrasekhar 
mass is reached before the BEG forms, so that collapse of the DM to a black hole occurs 
as soon as the BEG forms. The change in the slope of the curves around mx ~ 1 GeV is 
a combination of this effect with a decreased capture efficiency below mx ~ 1 GeV. We 
can see the formation of the BEG significantly improves the bound for light ADM. As an 
example, we see that for px = 1 GeV/cm^, BEG formation strengthens the constraint as long 
as mx ^ 13 GeV, while for higher masses the entire mass of ADM becomes self-gravitating 
before BEG formation occurs. 



C. Constraints from observed pulsars 



Now we consider the observations of a few relatively cold and old neutron stars that 
can provide tests of this effect. PSR J0437-4715 is a nearby pulsar at a distance of about 
139 ±3 pc from the solar system. The surface temperature is Tg = 1.2 x 10^ K |46| . When its 
secular motion is accounted for, calculations indicate its age is 6.69 x 10^ years |47|. Another 
nearby pulsar is PSR J2124-3358, located 270 pc away from us with a surface temperature 
Te < 4.6 X 10^ K [iel- Its age is 7.81 x 10^ years [4j 



For these nearby pulsars, we can calculate the central temperature from the surface 



temperature Te by using the analytical formula 48 

10"*^^ cm/s 



T ~ 1.288 X 10* K 



9s 



106 K 



0.455 



(42) 



where Qs = GMn/R^ is the surface gravity and we take ~ 1.7 x 10^^ cm/s^ for M„ = 
1.44 Mq and i?„ = 10.6 km. We find the central temperature is 2.1 x 10^ K and 2.5 x 
10^ K for J0437-4715 and J2124-3358, respectively. We take the ambient DM density to be 
0.3 GeV/cm^, because these pulsars are in our relative neighborhood. 

In Fig. 02]), we show the constraints on the DM-nucleon scattering cross section of the 
scalar ADM from the nearby pulsars J0437-4715 (left panel) and J2124-3358 (right panel). 
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Figure 2. Regions (colored) excluded by the nearby pulsars J0437-4715 (left) and J2124-3358 
(right). The shaded, diagonal and square cross-hatched, and black regions are as in Fig. ([T|). In 
hatched regions, the bounds are lifted. 



We can see that J0437-4715 can constrain scalar ADM with mx ^ 10 TeV without a BEC. 
With the formation of a BEC, it is also sensitive to the mass range rrix ~ 5 MeV — 13 GeV. 
The captured scalar ADM cannot form a BEC in the pulsar J2124-3358. This is because it 
has a relatively high central temperature, and the formation of a BEC requires a DM-nucleon 
cross section larger than the saturation cross section amax — 2.1 x 10~^^ cm^. 

Since the bound is sensitive to the DM density, we also consider neutron stars in regions 
with high px- Globular clusters possibly offer this type of environment, and observations 



of Pulsar B1620-26 place it in the globular cluster M4 49|] with an age of 2.82 x 10^ years 



46( 1 . Since it is far away from us, its surface temperature is unknown, and we are not able 
to calculate its central temperature. In our analysis, we take T = 10^ K as a reasonable 
approximation due to its advanced age. We take px = 10^ GeV/cm^ for the DM density and 



V = 20 km/s, motivated by discussions in Refs. 25|, |38|. Note that the exact value of DM 



density in globular clusters is unknown. Globular clusters are baryon-dominated systems 



and currently there is no evidence that DM is present in these systems; see Ref. [50[ for 
simulations of DM content in globular clusters. In Fig. ([3]), we show the constraints on the 
DM-nucleon scattering cross section of scalar ADM from the pulsar B 1620-26 in the globular 
cluster M4. Note that when the DM mass is larger than ~ 4.7 x 10^ GeV, Nbec ^ A^sei/ 
and all captured DM particles collapse before a BEC forms. 
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Figure 3. Regions (colored) excluded by the pulsar B1620-26 in the globular cluster M4. Note the 
globular cluster M4 is a baryon-dominated system, and there may in fact be no dark matter in this 
system. Here, we take px = 10^ GeV/cm^ motivated by numerical results in Refs. 



25 



. The 



shaded, diagonal and square cross hatched, and black regions are as in Fig. ([T]). In hatched regions, 
the bounds are lifted. 

VII. CONCLUSIONS 



We have studied the consequences of scalar ADM accumulation in neutron stars. Neutron 
stars have high density and are ideal objects for capturing DM at high rates. Since ADM 
does not self-annihilate, a high mass of DM can accrete in the neutron star, and, lacking 
Fermi degeneracy pressure, rapidly self-gravitate and exceed the Chandrasekhar limit. Fur- 
thermore, the formation of a BEC increases the density of the ADM by several orders of 
magnitude, which greatly accelerates the onset of gravitational collapse and considerably 
strengthens the constraints in certain regions of parameter space. For high dark matter 
mass, which corresponds to a low initial black hole mass. Hawking radiation can weaken the 
constraints. In the absence of light messenger particles, this effect is quite limited in scope, 
but if the black hole can heat the ADM directly it may affect a larger range of parameters. 

We have computed the size of all of these effects and found that some presently observed 
pulsars constrain scalar ADM far more tightly than what is currently possible with direct 
detection experiments. These constraints are stronger even than the upcoming generation 
of experiments. We also note that these constraints can be significantly improved in the 
future with observations of old pulsars in regions of DM density greater than 10^ GeV/cm^. 



Note added: Since the first version of our paper appeared, 5l| was written. In the 



case where a BEC forms and Hawking radiation is important, the authors of 5l| find that 
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scalar ADM with mass ~ 0.75 MeV — 16 GeV is excluded by nearby neutron stars if the 
scattering cross section with neutrons is sizable. Our results exclude scalar ADM in the 
mass range ~ 5 MeV — 13 GeV in the same scenario. The discrepancy in the low-mass limit 
comes about because we take into account the neutron degeneracy effects on the capture 
rate. This is important for DM mass below ~ 1 GeV. 
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